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Nonclassical rotational behavior at the vicinity of the A point 
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The rotational property of a quantum liquid at the vicinity of the A point T\ is examined. In a 
liquid helium 4 just above T\, under the strong influence of Bose statistics, the coherent many-body 
wave function grows to an intermediate size between a macroscopic and a microscopic one, which is 
of a different nature from the thermal fluctuations. It must reflect in the rotational properties such as 
the moment of inertia. Beginning with the bosons without the condensate, we make a perturbation 
calculation of its susceptibility with respect to the repulsive interaction, and examine how, with 
decreasing temperature, the growth of the coherent wave function gradually changes the rotational 
behavior of a liquid: The moment of inertia slightly decreases just above T\. This means that at the 
vicinity of T\, the mechanical superfluid density does not always agree with the thermodynamical 
one. We compare the result to the experiment by Hess and Fairbank. A new interpretation of the 
shear viscosity just above T\ is given from this viewpoint. 

PACS numbers: 67.40.-w, 67.40.Vs, 67.40.Db, 05.30.Jp 



I. INTRODUCTION 

A natural way to discuss superfluidity in a confined 
system is to focus on its rotational properties. When a 
liquid helium 4 is rotated at a temperature far above the 
A point T\, it makes a rigid-body rotation with a uniform 
vorticity rotv ^ owing to its viscosity. The angular 
momentum around z-axis has a form of Lf — I^fl, where 
if is a classical moment of inertia and il is a rotational 
velocity of a container. Figure. 1 schematically shows the 
f2 dependence of the angular momentum. When it is 
cooled to a certain temperature below T\, it shows two 
different behaviors according to the value of fi. Below 
the critical velocity il c , it abruptly stops rotating just as 
the system passes the A point, and rotv = is satisfied 
over the whole volume of the liquid. On the other hand, 
under a faster rotation than f2 c , the uniform vorticity 
abruptly concentrates to certain points, forming vortex 
lines and leaving other areas to satisfy rotv = 0. These 
phenomena instilled us with the notion that superfluidity 
abruptly appears at the A point. 

The basis of our phenomenological understanding is 
the two-fluid model, the foundation of which at T <C T\ 
has been established by microscopic theories. The basic 
assumption of the two-fluid model is that it completely 
separates the system into the normal and superfluid part 
from the beginning, and assumes that the latter abruptly 
emerges at T\ [l]. Furthermore, it assumes that the 
superfluid density defined in the mechanical properties 
completely agrees with that defined in thermodynamics. 

Here we must make a clear definition of superfluid- 
ity. Just above T\, anomalies in thermodynamical or 
mechanical constants of a liquid, such as the A shape 
of specific heat or the softening of sound propagation, 
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FIG. I: The schematic ^-dependence of the angular momen- 
tum L z . A solid straight line a is L z — ifQ, at T 2> T\, and 
a horizontal line b and a series of steps c is L Z (Q) at T < T\. 
A dotted curve d is a subject of this paper, a size of which 
is exaggerated for clarity. 



has been observed; These anomalies are normally consid- 
ered to arise from thermal fluctuations giving rise to the 
short-lived randomly oriented coherent wave functions. 
For superfluidity, however, the long-lived collective mo- 
tions of particles with a stable specific direction of motion 
are necessary. London stressed that superfluidity is not 
merely the absence of viscosity, but the occurrence of 
rotv = 0, and proposed an experiment to confirm this 
point J2| , which was later performed by Hess and Fair- 
bank 3] . This means that the complete disappearance of 
shear viscosity is attributed to rotv = over the whole 
volume of a liquid. Hence, the thermal fluctuations, be- 
ing a collection of independently growing and decaying 
wave functions, do not lead to superfluidity. (The relax- 
ation time of thermal fluctuations is far shorter than the 
characteristic time of the macroscopic rotational exper- 
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iments. Hence, they decay long before they affect the 
moment of inertia.) 

This paper will examine the foundation of the two- 
fluid model at the vicinity of the lamda point. Just 
above T\, particles experience the strong influence of 
Bose statistics, and therefore, the thermal-equilibrium 
coherent many-body wave function satisfying rotv = 
grows to a large but not yet macroscopic size. In contrast 
to thermal fluctuations, these long-lived wave functions 
have a possibility of affecting a mechanical property of 
the whole liquid, such as the moment of inertia. Specifi- 
cally, they may slightly reduce the moment of inertia I z 
just above T\, a possible L z (Sl) of which is schematically 
illustrated by a dotted line d in Fig.l. This is quali- 
tatively different from precursory phenomena owing to 
thermal fluctuations in that it requires a stable specific 
direction of motion. If it is true, it allows us to rede- 
fine the superfluid density in the mechanical phenomena. 
This means that at the vicinity of T\, the mechanical 
superfluid density does not always agree with the ther- 
modynamical superfluid density. 

The rotational properties of a liquid helium 4 has been 
subjected to considerable experimental and theoretical 
studies These studies, however, mainly focus on 

the dynamics of the quantized vortices in the superfluid 
phase in situations where the rotational velocity is not so 
small that the number of vortices are large. After the pi- 
oneering work by Hess and Fairbank Q , and by Packard 
and Sanders [f| , the regimes in which only a few vortices 
are present have rarely been explored. Hence, it is not 
surprising that almost no precise measurement has been 
made on I z just above the A point. 

A similar reason exists in theoretical studies as well. 
In theories of superfluidity, the infinite-volume limit is 
often assumed. In V — > oo, the only significant dis- 
tinction between states is that between the microscopic 
and macroscopic one, and therefore there is no room for 
the intermediate-sized wave functions in the theory. (In 
V — > oo, "large but not yet macroscopic" is substan- 
tially equivalent to "microscopic" p.) This clear-cut 
distinction lies behind the two-fluid model, and leads 
us to the preoccupied notion that superfluidity appear 
in a mathematically discontinuous manner at the Bosc- 
Einstein condensation temperature Tbec- For the real 
system, however, the overall transformation occurs more 
or less continuously. For the rotation of confined system, 
one cannot ignore the existence of the center of rotation 
and the boundary of the system, and therefore, one must 
take into account the size of the system. Hence, the va- 
lidity of the limit V — > oo is worth examination, and the 
magnitude of phenomena hidden in the V — > oo limit 
must be estimated by experiments. 

(1) When one views previous experiments of a liquid 
helium 4 from this point, one notes in the data by Hess 
and Fairbank an experimental sign suggesting a slight 
decrease of I z just above T\ (see Sec.4A). 

(2) The trapped atomic Bose gas opens a new possi- 
bility of precise measurements of I z just above T c . The 



measurement of the angular momentum using the pre- 
cession of a Bose-Einstein condensate of 87 Rb atoms was 
performed Q, which is analogous to the experiment by 
Hess and Fairbank, and by Packard and Sanders in a bulk 
liquid helium 4. It gives a data of L Z (Q) like b and c 
m Fig.l (Fig. 2 of Ref.0). Since the trapped Bose gas is 
a small system (= Ajim), its fl c is 10 4 times larger than 
that of a liquid helium 4 in Ref. [3| and [5| , which enables 
us to realize a situation in which only a few vortices are 
present. Whereas studies in this field are now centered on 
the anomalous behaviors at T -C T c , it has the potential 
for showing a slight change of I z just above T c . 

Although the deviation of the moment of inertia I z 
just above T\ from its classical value may be small, the 
essence of superfluidity is revealed in a primitive form in 
such a regime, which constitutes the necessary condition 
for discriminating a quantum fluid from a classical fluid. 
To consider these problems, we will make a somewhat 
different approach from conventional ones. At T < T\, 
the existence of the macroscopic coherent wave function 
4>{r) — \<f>\ exp[iS{r)\ (r is a center-of-mass coordinate 
of many helium 4 atoms) leads to rotv = geometri- 
cally, because the condensate momentum p is expressed 
by P = (H/myVS. Since we will focus on the continu- 
ous change of the system around T\, we cannot assume 
from the beginning the sudden emergence of <f>(r) at T\. 
Rather, considering the Bose system above and below T\ 
on a common ground will enable us to study the intricacy 
underlying the onset of superfluidity. We begin with the 
Bose system without the condensate, make a perturba- 
tion calculation of its susceptibility with respect to the 
repulsive interaction by taking peculiar graphs reflecting 
Bose statistics, and examine how the formation of the 
coherent wave function gradually changes the rotational 
behavior of the system. As a result, we derive a nonclas- 
sical rotational behavior that we normally think comes 
from rotv = 0, without assuming rotv = from the be- 
ginning. Specifically, we will derive the decrease of I z 
like d in Fig.l in the rotating repulsive Bose system just 
above T\ , and compare it with the data in a liquid helium 
4. 

This paper is organized as follows. Section 2 recapitu- 
lates the definition of the moment of inertia, and explains 
the physical reason of the nonclassical rotational behav- 
ior. Section 3 develops a formalism of the linear response 
of the system. (For the difference between this formalism 
and the thermal fluctuation theories, see Appendix. A.) 
Using the result in Sec. 3, Sec. 4 re-examines the exper- 
iment by Hess and Fairbank, and estimates the size of 
the intermediate-sized wave function and the strength 
of repulsive interaction in a liquid helium 4. Section 5 
considers the nonlinear response. From this viewpoint, 
Sec. 6 gives a new interpretation of the observed decrease 
of shear viscosity just above T\ in a liquid helium 4, and 
discusses some other examples. 
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FIG. 2: Schematic pictures of a part of the rotating bosons 
in a cylindrical container. White circles represent a initial 
distribution of particles. The rotation by long arrows moves 
white circles on a solid-line radius to black circles, (a) At 
T < Tx, the permutation symmetry holds over the whole 
liquid, (b) Just above Tx, it holds only within limited areas 
enclosed by a dotted lines, a size of which is exaggerated for 
clarity. 



II. MOMENT OF INERTIA OF THE 
REPULSIVE BOSE SYSTEM 

A. Moment of inertia 

Consider bosons in a uniform rotation around z-axis. 
For a liquid helium 4, the repulsive particle picture is not 
so unrealistic as it would be for any other liquid. Hence, 
as its simplest model, we use 

h=j2 e(p)& p $ p +uj2Yl K-«K>+ q ®p'*p> ( u > °)> 

V p,p' 1 

w 

where $ p denotes an annihilation operator of a spinless 
boson. 

The hamiltonian in a coordinate system rotating with 
a container is H — fi • L, where L is the total angular 
momentum. The rotation is equivalent to the application 
of a probe acting on a sample 8] [9]. The perturbation 
H ex = — fl ■ L is cast in the form — J2i(^ x r) • p , in 
which f2 x r = Vd(r) serves as the external field. Fig- 
ure. 2 shows a part of the Bose system in a cylindrical 
container. When the origin of r is put on the center 
of rotation, Vd(r) has a concentric-circle structure illus- 
trated by curved arrows in Fig. 2. (In the rigid-body ro- 
tation, Vd(r) agrees with the drift velocity at point r.) 
We define a mass-current density J(r), and express the 



perturbation H ex as 

-ft L = - J v d (r) ■ J(r)d 3 x. (2) 

Because of divvd(r) — 0, Eq.(2) says that Vd(r) acts 
as a trans verse- vector probe to the excitation of bosons. 
This fact allows us the formal analogy that the response 
of the system to Vd(r) is analogous to the response of the 
charged Bose system to the vector potential A{r) in the 
Coulomb gauge. Hence, J(r) in Eq.(2) has the following 
form in momentum space being similar to that in the 
charged Bose system 

J^r)=Y J (p+\) $l% +q e-^ T , (3) 

p,n 

(h = 1 and t = it). Vd(r) is a macroscopic external 
field causing the spatial inhomogeneity in the container, 
whereas J(r) contains both microscopic and macroscopic 
informations of the system. 

As the simplest susceptibility to Vd(r), we often use the 
mass density p = nm (n is the number density of par- 
ticles) as J(r) = pvd(r). Microscopically, however, one 
must begin with the generalized susceptibility consisting 
of the longitudinal and transverse part (p = x, y, z) 

X »v(q,u,) = 5^^( g , w )+ ^ X T {q,u). (4) 

By definition, the mass density p is a longitudinal re- 
sponse to an external force, p = x L (0i0)- As illus- 
trated in Fig. 2, however, the rotational motion of parti- 
cles is perpendicular to the radial direction, along which 
the influence of the wall motion extends into the con- 
tainer. Hence, in principle, one must use the trans- 
verse susceptibility x T (g,o;) for Vd(r) such as J(r) = 
[lim^o X T (<7,0)] v d (r). 

Using fl = (0, 0, O) in the left-hand side of Eq.(2) and 
the above J{r) and Vd = £1 X r = (— fly, fix, 0) in its 
right-hand side, one obtains the angular momentum L z 
as 

L z =x T (0,0) / (x 2 +y 2 )d 3 x-n. (5) 
Jv 

In a normal fluid, the susceptibility satisfies x T (0,0) = 
X^O, 0), and therefore the ordinary use of p is justified. 
The classical moment of inertia is given by 

lf=mn [ (x 2 +y 2 )d 3 x = X L (0,0) f (x 2 +y 2 )d 3 x. (6) 
Jv Jv 

In a superfluid, however, the above argument must be 
altered. For the later use, we define a term proportional 
to in by x^ 

/ \ x t, \ , fx L {q,u) - x T (q,u)\ 
Xnu{q,^) = <v„x (q, + q^Qu I ^ I 

= <Wx T (<7,w) +x liV {q,w), (7) 
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where x^v represents the balance between the longitudi- 
nal and transverse susceptibility [ioj . Comparing Eq.(5) 
to Eq.(6), one writes the moment of inertia I z — L z /fl 
using Xfiu 



i — r 



1 - - lim 

p 9^0 



-Xnu(q,o) 



(8) 



For the occurrence of nonclassical moment of iner- 
tia, the balance between the longitudinal and transverse 
low-energy excitation must be destroyed. In Eq.(8), 
lim q ^o[(q 2 / q ti q l ,)x fJ ,v(q, 0)] corresponds to the superfluid 
density p s . 

Consider x^f °f the ideal Bose system. Within the 
linear response, it is defined as 
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X$(?>0 = ^ J dTrapCiWnrXOl^J^feTj^Cff.OJlO), 

(9) 

where |0 > is the ground sate of ^2 e(p)& p & p . The term 
proportional to q^,q v has a form of 

LU + e{p) _ e{p + q) < ( 10 ) 

where f(e(p)) is the Bose distribution. 

(1) If bosons would form the condensate, f(e(p)) in 
Eq.(lO) is a macroscopic number for p = and nearly 
zero for p 0. Thus, in the sum over p in the right-hand 
side of Eq.(lO), only two terms corresponding to p = 
and p — —q remain, with a result that 



*#(<fc0) = Ps (T) 



(11) 



where p s {T) = mn c (T) is the thermodynamical super- 
fluid density, and n c (T) is the number density of parti- 
cles participating in the condensate. Equation. (8) with 
Eq.(ll) leads to 



h = r: i 



Ps(T) 



(12) 



(2) When bosons form no condensate, the sum over p in 
Eq.(lO) is carried out by replacing it with an integral, and 
one notes that q~ 2 dependence disappears in the result. 
Hence, using such a ic^v{<l,oj) m Eq.(8) leads to I z — 
if at q — > 0. This means that, without the interaction 
between particles, BEC is the necessary condition for the 
nonclassical moment of inertia. 

Under the repulsive interaction, however, the above 
argument is seriously affected. To see this, we must begin 
with a physical argument. 



B. Bose statistics and repulsive interaction 

We have a physical reason to expect the decrease of 
the moment of inertia in bosons at low temperature. 



The relationship between the low-energy excitations and 
Bose statistics dates back to Feynman's argument on the 
scarcity of the excitation in a liquid helium 4 [ill ], in 
which he explained how Bose statistics affects the many- 
body wave function in configuration space. To the rotat- 
ing bosons, we will apply his explanation. 

(1) In Fig. 2 (a), a liquid (white circles) is in the BEC 
phase, and the wave function has permutation symme- 
try everywhere in the container. Assume that the rota- 
tion of a container (depicted by curved arrows) moves 
white circles on a solid-line radius to black circles on a 
one-point-dotted- line radius (a transverse excitation) . At 
first sight, these displacements seem to be a large-scale 
configuration change, but this result is reproduced by a 
set of slight displacements (depicted by short thick ar- 
rows) from positions in the initial configuration to black 
circles after rotation. For any particle after rotation, it 
is possible to find a particle being close to it in the initial 
configuration. In Bose statistics, owing to permutation 
symmetry, one cannot distinguish between two types of 
particles after rotation, one moved from the neighboring 
position by the short arrow, and the other moved from 
distant initial positions by the long arrow. Even if the 
displacement made by the long arrows is a large displace- 
ment in classical statistics, it is only a slight displacement 
by the short arrows in Bose statistics. 

Let us imagine this situation in the 3N-dimensional 
configuration space. The above feature of Bose statistics 
means that in the configuration space, the excited state 
driven by rotation lies close to the ground state. Since the 
excited state is orthogonal to the ground state, the wave 
function corresponding to the excited state must spatially 
oscillate. Accordingly, the many-body wave function of 
the transversely excited state oscillates within a small 
distance in configuration space. Since the kinetic energy 
of the system is determined by the 3N-dimensional gra- 
dient of the wave function, this steep rise and fall of the 
amplitude means that the energy of the transverse exci- 
tation is not small even at q = 0, leading to the scarcity 
of the low-energy transverse excitation. This is the rea- 
son of x T (9jO) — * at q — ► below T\, whereas the 
particle conservation asserts that x L (it 0) = p is valid 
both above and below T\. Hence, X M1 y(g,0) in Eq.(7) 
changes to pq^q v /q 2 at q — ► 0, leading to I z — in 
Eq.(8). (This mechanism underlies the geometrical con- 
dition rotv s = 0.) 

(2) At high temperature, the coherent wave function 
has a microscopic size. If a long arrow of Vd(r) takes a 
particle to a position beyond the coherent wave function 
including that particle, one cannot regard the particle af- 
ter rotation as an equivalent of the initial one. The mech- 
anism below T\ does not work for the large displacement 
extending over two different wave functions. Hence, we 
obtain X T {q, 0) = x L {l, 0) at q ->■ 0, and I z = if. 

(3) Figure. 2(b) shows the boson system at the vicin- 
ity of T\ in the normal phase, in which the coherent 
many-body wave function grows to a large but not yet 
a macroscopic size (regions enclosed by a dotted line). 
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The permutation symmetry holds only within each of 
these regions. When particles are moved from a region A 
to another region A', the mechanism below T\ does not 
work. 

The repulsive interaction U between particles, how- 
ever, affects this situation. In general, when one moves a 
particle in the interacting system, it induces the motions 
of other particles. In particular, the large- distance dis- 
placement of a particle in coordinate space causes the 
excitation of many particles, and therefore it needs a 
large excitation energy. This means, in the low-energy 
excitation of the system, one observes mainly the short- 
distance displacement of particles. When applying this 
tendency to the low-energy excitation of repulsive bosons, 
one knows that excited particles are not likely to go be- 
yond a single coherent wave function, but likely to remain 
in it, and therefore the mechanism working below T\ 
works just above T\ as well. This view will be tested as 
follows. If we increase the strength of U in x M „(g, 0), the 
excited bosons get to remain in the same coherent wave 
function, and therefore the low-energy transverse excita- 
tion will raises its energy owing to Bose statistics as dis- 
cussed in (1). Hence, the condition of x L (li 0) = X T (<7> 0) 
at q — > will be violated at a certain critical value of U. 
Alternatively, if we decrease the temperature at a given 
U, the above condition will be violated at a certain tem- 
perature T on . 

A geometric feature inherent in the rotation will play 
an important role in the above mechanism. The external 
field Vd{r) has the structure of concentric circle; hence, 
the center of rotation is a fixed point. The displacements 
of particles near the center is so small that they do not 
go beyond a single coherent wave function (a region B 
in Fig. 2(b)). The center of rotation is the most probable 
point for the mechanism discussed in (a) to work. Hence, 
the region near the center is most likely to decouple from 
the motion of container. With decreasing temperature, 
this decoupling will extend from the center to the wall, 
which depends on the rotational velocity 0. 

To formulate these mechanisms, we consider the per- 
turbation expansion of \tiv with respect to U in Sec. 3. 
After comparing it with the experiment in Sec. 4, we ex- 
tend it to the nonlinear response to Vd(r) in Sec. 5. 



III. LINEAR RESPONSE 

We will formulate the moment of inertia in the repul- 
sive Bose system at the vicinity of T\. In the integrand 
of Eq.(9), one must use, instead of |0 >, the ground state 
\G> of Eq.(l) as follows 



(G\T t J^x,t)J„(0,0)\G) 



(0\T T J tl (x,T)X(0,0)exp 


\~f 

Jo 


dr Hi (t ) 


10) 

i 


(0\exp 


Jo 


10} 



(13) 




p+q 



p+q 



(0 



FIG. 3: (a) The first-order Feynman diagram of a current- 
current response tensor J M J„ (a lower bubble), and an exci- 
tation owing to the repulsive interaction U (an upper bubble 
with a dotted line). The black and white small circle repre- 
sents a vector and scalar vertex, respectively. The exchange 
of particle lines between the excitation and the response ten- 
sor yields (b). Similarly, the interchange of particle lines in a 
deformed square yields (c). 



where Hj(t) represents the repulsive interaction. Fig- 
ure 3(a) illustrates the current-current response tensor 
Jfj,(x, t) J„(0, 0) (a lower bubble) in the medium: Owing 
to exp(~ f Hi(t)cIt) in Eq.(13), scatterings of particles 
frequently occur in \G > as illustrated by an upper bub- 
ble with a dotted line U in Fig. 3 (a). (The black and 
white circle represents a coupling to Vd(q) and to U, re- 
spectively.) The state \G > includes many interaction 
bubbles like the upper one in Fig. 3 (a) with various mo- 
mentums p' and p' + q'. A solid line with an arrow rep- 
resents 



G(iuj n ,p) 



1 



e(p) 



(14) 



(/i is a chemical potential implicitly determined by 
V-^exp^ete) + £ - n]) - l]" 1 = n). Owing to the 
repulsive interaction, the boson has a self energy E (> 0) 
(we ignore its u> and p dependence by assuming it small). 
With decreasing temperature, the negative \i at high tem- 
perature approaches a small positive value of E, finally 
reaching Bose-Einstein condensation satisfying fi = E. 

When the system is just above T\ in the normal phase, 
particles in the ground state \G > are under the strong in- 
fluence of Bose statistics. Hence, the perturbation must 
be developed in such a way that, as the order of the 
perturbation increases, the susceptibility gradually in- 
cludes a new effect owing to Bose statistics. Specifically, 
the lower bubble J M (x, r) J„(0, 0) and the upper bubble 
in Fig. 3(a) form a coherent wave function as a whole. 
When one of the two particles in both the lower and 
upper bubble have the same momentum (p — p'), and 
the other in both bubbles have another same momentum 
(p + q = p' + q') in Fig. 3(a), a graph made by exchang- 
ing these particles must be included in the expansion. 
Such a transformation in Fig. 3(a) takes place as follows. 
The exchange of two particles having p and p'(= p) by 
thick white arrows yields Fig. 3(b). Furthermore, the in- 
terchange of two particles having p + q and p + q' {= p + q) 
in Fig. 3(b) by thick white arrows yields Fig. 3(c). The re- 
sult is that two bubbles with the same momentum arc 
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linked by the repulsive interaction, the contribution of 
which to Xfiv is given by 



4i> 



Is I , 



/(e(p) + E)-/(e(p + <?)+£) 



w + e(p) - e(p + q) 



(15) 



(a) With decreasing temperature, the coherent wave 
function grows to a large size, and the interchange of 
particles owing to Bose statistics like Fig. 3 occurs many 
times. Hence, one cannot ignore the higher-order terms 
in Eq.(13), which become more significant with the 
growth of the coherent wave function. 

(b) Among many particles contributing to Eq.(15), 
particles stationary to a container play a dominant role. 
Specifically, a term with p = in Eq.(15) corresponds to 
an excitation from the rest particle, and that withp = — q 
corresponds to a decay into the rest one. 

These two considerations (a) and (b) lead to the fol- 
lowing form of xj}} (q, 0) at the vicinity of T\ 



xLl } M) 



q^q v 



2 v 



(16) 



1=0 



where Fp(q) is 

(exp(/3[E - a*]) - I)" 1 - (exp(/J[e(g) +£-/*])- 1)" 



(17) 



a positive monotonously decreasing function of q 2 , which 
approaches zero as q 2 — > oo. 

At a high temperature (f3fi -C 0) in which Fp(q) is 
small, a small Fp(q) guarantees the convergence of an 

infinite series in xijv(c[, 0) °f Eq.(16), with a result that 



2 Vl-UFp(q) 



(18) 



With decreasing temperature, however, the negative /x 
gradually approaches E, hence £ — \i — > 0. Since Fp(q) 
increases as £ — /j — * 0, it makes the higher-order term 
significant in Eq.(16). An expansion form of Fp(q) = 
F/3 (0) — aq 2 + ■ ■ ■ around q 2 = has a form such as 



ft 



2 /|/3[ M (T)-£] 



4sinh 



1 



1_£ 

2 . , /|^(T)-E]|\ 2m 



tanh 



.(19) 



In E — /i — > 0, UFp{0) increases and finally reaches 1, 
that is, 



U(3 = 4sinh 



2 (0MT)--E(U)] 



(20) 



At this point, the denominator in the right-hand side 
of Eq.(18) gets to begin with q 2 , and x < j}}(q,0) therefore 
changes to a form of q^q v I q 2 at q — > 0. This means that a 
non-zero coefficient Fp{Q)/{2VUa) of q^q^/q 2 appearing 
in Eq.(18) gives a non-zero value of x L (<iS 0) — x T (<Zi 0) in 
Eq.(7), hence the moment of inertia shows the nonclas- 
sical behavior in Eq.(8). 

From now, we call T satisfying Eq.(20) the onset tem- 
perature of the nonclassical moment of inertia T on [l^] • 

At the vicinity of T\, Eq.(20) is approximated as Uf3 = 
/3 2 [/x(T)-E([/)] 2 for a small /x-E. This condition has two 
solutions n(T) — E([/)±v / C^bT. It is generally assumed 
that the repulsive Bose system undergoes BEC as well as 
a free Bose gas. Hence, with decreasing temperature, 
fx(T) in the presence of repulsive interaction U should 
reach E(f7) at a finite temperature, during which course 
the system necessarily passes a state satisfying fJ,(T) = 
£(£/) — y/V ksT [13J. One concludes that the nonclassical 
rotational behavior always occurs prior to BEC in the 
repulsive bosons , that is, T on > T\. 

The chemical potential /1, hence /i — E as well, deter- 
mines the size of the coherent many-body wave function 
[lij (l5j , which corresponds to the size of regions enclosed 
by a dotted line in Fig. 2(b). The emergence of q~ 2 singu- 
larity in xfiv (<?, 0) in the process of E — fj, — ► is a math- 
ematical expression of the instability mechanism induced 
by the growth of the coherent wave function. When 
U is small, this instability occurs after the wave func- 
tion grows to a large size corresponding to a small E — ji. 
When U is large, this instability already occurs at a larger 
E — [i in which the wave function is smaller than the for- 
mer one. 

At the onset temperature T on , substituting Eq.(19) 
into Eq.(18), we find x$ at q — > 

*£>(?, 0) = #4 tanh ( i/WM(T r )-E]h ^ 



and with the aid of Eq.(20) 



q^qv 



^sinh|^ „[/i(r on )-E]| q 2 



X$ (?> 0) is given by 



<?m q v 



X$(Q>°) = mc (T<m)no(T on ) ^ 2 



(21) 



(22) 



(23) 



where 



At q — > 0, the denominator 1 — UFp(q) in the right- 
hand side of Eq.(18) has a form of [1 - UFp(0)] + Uaq 2 . 



n (T) 



Vexp(-f3[fi(T)-E])-r 



(24) 
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is the number density of p = bosons, and 

2 



c(T) 



exp(/3HT)-£|) + l' 



(25) 



is a Fermi-distribution- like coefficient. For the finite 
system just above T\, n (T) has a large but not yet 
macroscopic value. In the theoretical limit V — > oo, this 
quantity is normally regarded to be zero. In real finite 
system, however, its magnitude must be estimated by ex- 
periments (see Sec. 4). Using Eq.(23) in Eq.(8), we obtain 



/:i7;,„) = /:' ll- c (T on )?^^W? 



1 - 



Ps(T : 



(26) 

where p s {T) = mc(T)no(T) is the mechanical superfluid 
density. 

p s (T) reflects the intermediate-sized coherent wave 
function. "Intermediate" means that it does not play 
the role of order parameter characterizing the thermody- 
namical phase, but affects mechanical properties of the 
system. At T = T on , the system shows a small but finite 
jump from if to I z , which is proportional to no(T on ). 
Since c(T on ) has an order of one, the magnitude of 
this jump is determined mainly by the non-macroscopic 
fk>(T on ) an d the moment of inertia therefore only slightly 
decreases from the classical value. (This does not mean 
that the thermodynamical quantities show a finite jump.) 
At T < T on , the moment of inertia varies with T follow- 
ing I Z (T) in Eq.(26). Equation (26) determines an initial 
slope of the dotted curve d at = in Fig.l. When the 
system reaches T = T x , the condition of p = £ makes 
p s (T x ) = Ps {T x ) because of c(T A ) = 1, n (T x ) = n c , 
which shows a natural connection of Eq.(26) to Eq.(12) 
with the thermodynamical p s (T). While the thermody- 
namical p s (T) satisfy p s = (m 2 k B T /h 2 )\<j)\ 2 at T < T x , 
one can expect no simple relation of the mechanical p s (T) 
with the thermodynamical quantities at T\ < T < T on . 



IV. COMPARISON TO EXPERIMENTS 

A. Experiment by Hess and Fairbank revisited 

Hess and Fairbank, after they confirmed that a su- 
perfluid 4 below T x remained at rest under the ex- 
tremely slow rotation, made another type of experiment 
in their classic paper Ref.[3[. First, at an initial tem- 
perature T\ below or above T\, they rotated a liquid 
helium 4, contained in a small cylinder of radius R = 
0.44 mm, at fl =1.13 rad/s. Later they heated it up 
to the temperature as it comes into rigid-body rotation, 
and precisely measured a small change A51 of the angu- 
lar velocity. Using Eq.(12), the rotational energy before 
heating is 1/2 x lf(l — p s /p)Q 2 , whereas it changes to 
1/2 x lf(Vl — Aft) 2 after heating. By the conservation 
of energy, the former must be equal to the latter, with 
a result that 2Af2/Q = p s /p. Below Tx, a fraction of 



liquid does not participate in the rotation, whereas af- 
ter heating it does. The rotation after heating therefore 
always becomes slower, and 2Af2/f2 = p s /p is positive. 
They plotted AQ /f2 as a function of initial temperature 
T\ (Fig.2 of Ref.[|). 

When the initial temperature Ti was lower than Tx, 
AVt I Q was properly explained by the theoretical value of 
p s (Ti)/p. If the two-fluid model is exactly valid near 
Tx, p s must vanish at T > T\. Hence, AQ/Q mea- 
sured under the condition of T\ > Tx must be exactly 
zero. For T x = T x + 0.03K and T x + 0.28K, however, 
the measured Afi/fi in Ref [3j were not exactly zero. 
Although the error bars were large compared with its 
absolute values, its central values were significantly dif- 
ferent from zero: AQ/D, = 4x 1(T 5 at T x = T x + 0.03K, 
and 1.5 x 10~ 5 at Tx + 0.28K. A natural interpreta- 
tion of this result is that these Ti's were lower than 
the onset temperature T on of the nonclassical moment 
of inertia, and at such Ti's, I z was already slightly 
smaller than its classical value. If it is true, instead of 
Eq.(12),7,(T) = lf(l-p s {T)/p) (Eq. (26)) must be used, 
and we obtain 2A£7/f2 = p s (T±)/p, Hence, we obtain 
p s {Tx + 0.03K)/p = 8 x 10- 5 , and p s [T x + 0.28K)/p = 
3 x 10 -5 Since the number density of atoms in a 

liquid helium 4 is n = 2.2 x 10 22 atoms/cm 3 , we have 
no = 10 18 atoms/cm 3 . This is an experimental estima- 
tion of the number of helium 4 atoms participating in 
the intermediate-sized coherent wave function in a bulk 
helium 4 just above T x . Two temperatures, Ta+0.03K 
and T>+0.28K, are situated within the temperature re- 
gion (Tx < T < 2.8 K) in which the viscosity begins to 
decrease above T x . In Ref.0], the author's focus was on 
the rigidity of a superfluid against the rotation (the re- 
sult was later named Hess- Fairbank effect), and they did 
not mention the small non-zero value of Ail/Q above T\. 
Although this point did not attract the interest of many 
people, it is worth studying closely in the future. 



B. Estimation of the repulsive interaction U 

Let us make an estimation of the repulsive interac- 
tion U using the approximate form of Eq. (20) , U (3 on — 
P 2 n (p — X) 2 . (1) For the present, we suppose T on = 2.8K 
by assuming that the decrease of viscosity just above Tx 
has a relation to the emergence of the intermediate-sized 
coherent wave function (see Sec. 6. A). (2) We assume that 
p(T) - £(Z7) in Eq.(20) follows the formula 



M (T) - E(J7) 



/ ff3/2(l) 

V 2v^ 



3/2 



(27) 

(g a (x) = J2 n a -" / n<1 ) 011 t ne assumption that the particle 
interaction U and the particle density of a liquid helium 
4 are renormalizcd to T x = 2.17 K (an approximation 
that dates back to London). Hence, we obtain a rough 
estimation of U as 0.5 x 10~ 17 erg. This value is approx- 
imately close to the repulsive interaction U c obtained by 
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Bogoliubov's spectrum c = (h/m)y^ma, but somewhat 
smaller than it. (The velocity of ordinary sound c = 220 
m/ s of a liquid helium 4 near T\ gives a scattering length 
a = 0.7 ran, hence {J c = h 2 /(ma 2 ) = 3.4 x 10~ 17 erg.) It 
is a difficult problem to relate these values to the realistic 
potential between helium 4 atoms such as the Lennard- 
Jones potential U(r) = 4e[(2.556/r) 12 - (2.556/r) 6 ] with 
e = 1.41 x 10" 15 erg. 



V. NONLINEAR RESPONSE 

When the precise measurement of I z is performed, the 
dynamic response of I z (fl) will become a next subject, 
which appears in the nonlinear response of the system. 
As discussed in Sec. 2, J contains both microscopic and 
macroscopic information, whereas Vd{r) is a macroscopic 
external field, and therefore the susceptibility connect- 
ing J and Vd appears as J(r) = [lim 9 ^ X T (<L 0)] v d{r). 
Hence, we consider not a general form of the nonlinear 
susceptibility, but a correction term x T ' non ( v d) to the lin- 
ear response, such as J = [x T '^ 1 ' + X T ' n ° n ( v d)]vd- Since 
X T ' non (vd) does not depend on the direction of Vd, it in- 
cludes only even powers of Vd = fir, which leads to the 
spatial inhomogeneity and the dynamic response. 

For the dynamic response, we define some quantities. 
The current J(r) = x T (0,0)vd(r) in Sec. 2. A is replaced 



by J(r) = x T (0, 0, r)v ( j(T"), where r—y/ x 2 + y 2 is a 
distance from the center of rotation. Correspondingly, 
instead of Eq.(5) and (8), we define 



and 
4 (ft) 



lim 

q^0 



X (0,0,fi,r)r¥i-n, 



-X^(q,0,n,r) 



%1v 



(28) 



r l d 6 x. (29) 



The position-dependent rotational velocity is defined as 



Oo(r) 



1 lim 

p 9^0 



<3V °» 



-x^(q,0,ft,r) 



ft. (30) 



We extend the mechanical superfluid density p s (T) in 
Eq.(26) so that it has ft and r dependence and satisfies 
p s (T,fl) = J p s (T,fl,r)2nrdrdz as 



p s (T,fl,r) = lim 



qp.q v 



-x^(q,0,ft,r) 



(31) 



Let us consider the first approximation of the above 
quantities. We begin with 



< J^x,t) >=< GISU^x^SIG >, 



where S 



Texp 



-i f^dt'Hexir,?) 



(32) 
Using 



H ex (r) = — v^(r)Jfj,(r), the analytical continuation t — > 
r = it is performed in the higher-order expansion terms 




(b) 




■o 



+q pl+q 



(c) pQ- 




FIG. 4: (a) The third-order Feynman diagram of the current- 
current response tensor accompanied by the bubble excita- 
tion. (The definition of symbols is same as in Fig.3.) The ex- 
change of particle at one vertex yields (b). Similar exchanges 
at the other vertices yield (c). 



in the right-hand side of Eq.(32) [l7|. As the simplest 
nonlinear susceptibility for J^, we consider the third- 

(3) 

order term Xp'v^rVdVd v d with respect to H ex and extract 
a correction term to the linear susceptibility x)u> ( v d)vd 

from X^cr^^d w d • X^( v d) is illustrated as a square 
in Fig. 4(a), in which we choose only one vertex out of 
three vertices for the coupling to Vd(q) and make the 
others couple to the coarse-grained external field Vd(r). 
The lowest-order dynamic response of I z (ft) is obtained 

(3) (3) 

through xW( v d)- (From now, we express xW( v d) by 
X$ 

Comparing Fig.3 (a) and 4(a), we find that p in the 
lower bubble of Fig. 3(a) splits into p and p\ in Fig. 4(a). 
We obtain a formula corresponding to the square in 
Fig.4(a) 



X 



n,m p,pi 

pi\ f p + pi 



2 / V 2 

x G(iu) n + iuj,p + q)G(iuj ni p) 

x G(iuj m + iu),pi + q)G(ico m ,pi), (33) 

where [(p + Pi)/2] • [(p + Pi)/2 + q] comes from the 
coupling to the upper and lower Vd{r) in Fig. 4(a). In 
general, a loop with four vertices has three inner frequen- 
cies u) n ,u> m , and u>i. For the susceptibility like Fig. 4(a), 
however, the q,u> in Xui/(q^^) enters at one of the four 
vertices and leaves at another, thus leaving only two fre- 
quencies u> n , uj m as internal ones in Eq.(33). The re- 
maining H ex in Jq drH ex , corresponding to u>i , appears 
as /3H ex (= /3no\vd(r)\ 2 ), since the macroscopic Vd(r) in 
H ex slowly varies with r. 
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(a) With decreasing temperature, the coherent wave 
functions gradually grow, and therefore the particle in- 
terchange owing to Bose statistics frequently occurs in 
Xfiu • By applying Eq.(13)-likc formula for the product 
of four currents, we obtain a perturbation expansion of 

(3) 

Xfiu with respect to the repulsive interaction Hj. Fig- 
ure. 4(a) illustrates the square surrounded by a bubble in 
\G > (an analogue of Fig. 3(a)). The particle interchange 
at one of three vertices yields Fig. 4(b). Fig. 4(c) shows a 
horizontal and vertical extension of bubble chains. 

(b) With decreasing temperature, particles stationary 
to a container get to play a dominant role. Compar- 
ing Fig. 3(c) and Fig. 4(b), we apply the argument above 
Eq.(15) to Eq.(33), knowing that p = and p — —q plays 
a dominant role in the sum over p. In the sum over pi 
as well, the dominant process comes from p\ = and 
Pi = —q. At first sight, there seems to be four combina- 
tions in (p,pi). Owing to (p+pi)/2 and (p+p 1 )/2 + q'm 
Eq.(33), however, only (p,pi) = (0,—q) and (— g,0) are 
possible. In Fig. 4(a), the coupling to Vd(q) is possible on 
the upper or lower vertex as well. Although this case has 

a different expression of x\iv (q, iw) from Eq.(33), the par- 
ticle interchange and the dominance of p = or p = — q 

particles derive the same formula from its xW \1i §i y " 
ing the final result the symmetry factor 2. Hence, one 
obtains the term representing the balance between the 
longitudinal and transverse susceptibility for the square 
in Fig. 4(a) 



\F p {q)\ 2 {nrfPn a (T). 

(34) 



The process from Fig. 4(a) to (c) changes Eq.(34) as 
follows. Using Eq.(13)-like formula for the product of 
four currents, bubble chains in Fig. 4(c) are extended to 
infinity. Among these terms, to derive the nonclassical 
I z (0) in Eq.(29), we pick up only terms that give a factor 
proportional to q^q v /q 2 at q — ► limit because other 
terms vanish in the q — > limit of Eq.(29). Hence, 



x$(<z,o,n,r) = - 



q^qv 



l*M<?)l S 



[l-UFp(q)Y 



2 

(nr) 2 f3n (T). 



1 1 



(35) 



Similarly to Eq.(18), the instability occurs in Eq.(35) 
when the condition of Eq.(20) is satisfied, because the 
denominator in the right-hand side of Eq.(35) gives q~ 4 , 
and therefore the coefficient of q^q v diverges as q~ 2 at 
q = 0. At T — T on , the same procedure as that from 
Eq.(18) to (21) yields 

1 / 2m \ 2 1 

X^(9.0,n,r) = ( ) l3on^n Q (T on ) 

x tanh 2 ( PonlKTon)-^ (fir) 2M 1) (36) 



and with the aid of Eq.(20) 

V A V k B T or . 



1 



1 



?m q v 



n 2 sinh 2 (/3 on [ M (T on )-S]) q 2 
Using p s (T) = mc(T)n a {T) as in Eq.(23), we obtain 



(37) 



Xm/ ) (?.°»°. r ) = -mn (T on ) 



Ps(T on )\ ( m(Q,r) 2 \ q^ 



kBT n j q 2 
(38) 



Substituting Eqs.(23) and (38) to Eq.(30) and using 
mno/p = c~ 1 {p s /p) from the definition of p s , we get at 
T<T on 



n (r) 



Ps(T) 



P 



=(T) 



p s {T)Ym(nr) 2 



P 



k R T 



n.(39) 



Using Eq.(39), we obtain a velocity field vo(r) — tto(r) x 
r, and a vorticity field 



rotvo(r) = 2fl 



Ps(T) 
P 



c(T) 



Ps(T)Y m(Or) s 



P 



k R T 



(40) 

This form shows a little change from rotv — 2fle z to 
rotv = 0, that is, from a normal fluid to a superfluid. 
These results are summarized in the mechanical super- 
fluid density 



p.(T,n,r) = p s (T) 



1 



=(T) 



p s (T)\ 2 m(n r y 



k B T 



10- 



(41) 
sec 2 /cm 2 . 



In a liquid helium 4, m/ks — 2.13 x 
Hence, the dependence of Eq.(41) is negligibly small. 

When the container is rotated in the normal phase, a 
liquid makes the rigid-body rotation satisfying Vd(r) = 
fi x r. In the superfluid phase, a quantum vortex appears 
and a liquid makes a peculiar rotation satisfying v(r) = 
{hl/m)e z x r/r 2 (I is an integer) being independent of 
fL At T\ < T < T on , however, while a liquid follows 
the rotation of a container, a region around the center of 
rotation slightly reduces its rotational velocity. Hence, 
the velocity field in a container takes an intermediate 
form between the velocity near the center and at the 
boundary. One can call it a differential rotation (l8| . 

(1) Near the boundary, Vd(r) = 1~2 x r becomes a large 
external field, and more higher-order terms than x^ will 
contribute to I Z (Q). Presumably, this will suppress su- 
perfluidity, and vq(t) = Oo(r) x r with Eq.(39) will ap- 
proach Vd{r) = Jl x r of the rigid-body rotation. 

(2) When one increases the rotational velocity of a 
container, one can expect a similar result. With increas- 
ing Q, the region rotating more slowly than the container 
will shrink to the center of rotation, and therefore the 
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FIG. 5: The temperature dependence of the shear viscosity 
of a li quid helium 4. (The data at T > T\ is taken from Fig. 3 
of Ref.[l9j. The data at T < T\ that strongly depends on the 
experimental method is written schematically.). 



ent wave function must affect the response of the system. 
One must have a microscopic theory describing the in- 
fluence of superfluidity on the coefficient of viscosity r], 
irrespective of whether the decrease of rj is caused by the 
fluctuations or by the long-lived wave functions. For such 
a theory (a) one must apply the linear-response theory 
not to the mechanical, but to the thermal perturbation 
[20l | . The formulation of the latter perturbation includes 
more subtle points than that of the former one owing to 
the thermal dissipation, (b) The mechanism of shear vis- 
cosity in a liquid, which has similarities with the motion 
of dislocations in a solid, is considerably different from 
that of shear viscosity in a gas. Hence, the quasiparticle 
approximation, which is often used for the excitations in 
a liquid helium 4, is questionable for the shear viscosity 
in a liquid. The microscopic theory of the influence of 
superfluidity on 77 is a future problem. 



nonclassical I z will approach the classical one as illus- 
trated by a dotted line d in Fig.l. 

(3) As T — > T\, the region near the center of rota- 
tion, which has a smaller vorticity, will enlarge toward 
the wall. (To describe such a growth and shrink, the ef- 
fect of the wall must be taken into account as a boundary 
condition.) 

(4) At T < Ta, the nonclassical behavior discussed in 
this paper is masked by the quantum jump of L z owing 
to the emergence of the quantum vortex. If the quantum 
vortex could be suppressed, one would see the hidden 
behavior of L z (fi), that is, Eq.(39) with p s (T)/p = 1. 



VI. DISCUSSION 
A. Shear viscosity at the vicinity of T\ 

In a classical liquid, with decreasing temperature, the 
shear viscosity gradually increases. The shear viscosity of 
a liquid is inversely proportional to the rate of a process 
in which a hole in a liquid propagates from one point 
to another over the energy barriers. With decreasing 
temperature, this rate decreases, thus increasing the vis- 
cosity. In a liquid helium 4, however, it reaches a maxi- 
mum value at 2.8K, and begins to reduce its value, finally 
dropping at the A point (2.17K). Figure. 5 illustrates this 
temperature dependence taken from the data at T > T\ 
111 LL9JI. In view of some precursory phenomena observed 
in a liquid helium 4, it is natural to explain it in terms 
of the thermal fluctuations giving rise to the short-lived 
wave function. 

From the viewpoint of this paper, however, one must 
point out two other aspects of this phenomenon. (1) 
Under the strong influence of Bose statistics just above 
T\, the large but not yet macroscopic long-lived coher- 



(2) The viscosity appears in the Navier-Stokes equation 
as r/Av, which has a form of —rj ■ rot(rotv) in an incom- 
pressible fluid. The decrease of the viscosity comes from 
either the decrease of 77 or that of rotv 0] ■ The former is 
owing to the microscopic change of the system, whereas 
the latter arises from the macroscopic transformation of 
the velocity field v(r) in a container. At T < T\. the 
temperature dependence of the shear viscosity strongly 
depends on the experimental methods, such as (1) a 
Poiseuille flow method (a lower horizontal line at T < T\ 
in Fig. 5), (2) an oscillating disc viscometer method (an 
upper curve at T < T\) and (3) a rotation viscometer 
method. Since the microscopic mechanism determining 
77 is not sensitive to the types of the macroscopic mea- 
surements, this fact suggests that the shear viscosity at 
T < T\ is strongly dependent on the macroscopic change 
of rotv induced by the measurements. Above T\, how- 
ever, since the coherent wave function has not yet grown 
to a macroscopic size, it is natural to attribute the de- 
crease of viscosity at first to that of 77. Just above T\, 
however, we have another possibility in —r)-rot(rotv) that 
the emergence of the intermediate-sized coherent wave 
functions contributes to the decrease of rotv, because 
they are irrotational. This means that the microscopic 
mechanism of a liquid affects not only the coefficient 77, 
but also the macroscopic velocity field v(r). This is a 
peculiar case in fluid mechanics. In fluid mechanics, the 
velocity field v(r) is a solution of the equations of motion 
with given coefficients made of mechanical or thermody- 
namical constants of a liquid. In the system such as a 
liquid helium 4, one cannot clearly distinguish between 
the microscopic and macroscopic phenomena, and there- 
fore, this clear separation of coefficient and solution of 
the equation is not obvious in contrast with a classical 
liquid. One must make a different approach that will ex- 
amine the foundation of the two-fluid mechanics at the 
vicinity of T\. 
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B. Various manifestations of superfluidity 

Superfluidity is a complex of phenomena, and there- 
fore has some different manifestations, such as (1) per- 
sistent current without friction, (2) the Hess-Fairbank 
effect, (3) quantized circulation, (4) almost no friction 
on moving objects in the system below the critical veloc- 
ity, (5) peculiar collective excitations and (6) the Joseph- 
son effect. The conventional thermodynamical definition 
of superfluid density p s {T) has been proved to be use- 
ful for describing various mechanical manifestations of 
superfluidity. The result of this paper implies that the 
superfluid density in the mechanical phenomena does not 
always agree with the thermodynamical p s (T), and that 
the interplay between Bose statistics and the repulsive 
interaction sometimes require us to define the mechani- 
cal superfluid density p s (T) . This p s (T) may have some 
different definitions, each of which is specific to each man- 
ifestation of superfluidity, such as p s (T) in Eq.(26) for the 
rotation. Hence, we must consider the existence or non- 
existence of p s (T) in each manifestation on a case-by-case 
basis. Here we make a comment on three examples. 

(A) The Meissner effect in the charged Bose system is a 
counterpart of the nonclassical rotational behavior in the 
neutral Bose system [8j [9( . (The enhanced diamagnetism 
above Tg ec appears as the precursory phenomenon ow- 
ing to thermal fluctuations |2l|.) The phenomenon like 
the decrease of I z , not owing to fluctuations but requiring 
a transformation extending to the whole system, occurs 
just above Tbec m the Meissner effect as well: With de- 
creasing temperature, the charged Bose gas with short- 
range repulsion begins to exclude the magnetic field prior 
to the BEC [22j. In contrast with the enhanced diamag- 
netism, the exclusion of the applied magnetic field just 
above Tbec implies that the large but not yet macro- 
scopic wave functions cause a change extending to the 
whole system. This means that Bose statistics is essential 
for the Meissner effect, but the existence of the macro- 
scopic condensate is not the necessary condition. The 
same form of the mechanical superfluid density p s {T) is 
useful for the Meissner effect as well. 

(B) The nonclassical rotational behavior above T\ will 
become more realistic in superfluidity of small systems. 
Recently, a helium 4 droplet consisting of about 10 4 he- 
lium 4 atoms is found to show a sign of superfluidity. 
The infrared rotational spectrum of small molecules at- 
tached to the helium 4 droplet, such as oxygen carbon 
sulfide, shows a signal indicating a significant change of 
its moment of inertia around a molecular axis at T < T\, 
which suggests a transition occurring in the surrounding 
helium 4 environment [23j. Although the experimental 
condition, such as the temperature or the rotational ve- 
locity, cannot be freely controlled until now as in the bulk 
helium 4, it will open a new probability of superfluidity of 
small systems. If the same experiment as in Sec. 4 could 
be made for a liquid helium 4 droplet, one would see a 
larger role of the intermediate-sized coherent wave func- 
tion p s (Ti)/p in the rotational properties @. When the 



number of atoms in a droplet is too small, however, the A 
transition will become obscure. Hence, there may be an 
optimum size of the system for detecting the nonclassical 
rotational behavior just above T\. 

(C) Solid helium 4 has been termed a quantum crys- 
tal. Recently, an abrupt drop in the moment of inertia 
was found in the tortional oscillation measurements on 
solid helium 4 confined in a porous media [24| and on a 
bulk solid helium 4 [25[ . This discovery leads us to recon- 
sider the definition of superfluidity and that of solids [2(| . 
The fundamental feature of crystals is their periodicity in 
density; that is, diagonal long-range order (DLRO). One 
has to face a serious question whether crystals remain 
stable while showing superfluidity that violates their pe- 
riodicity. Phcnomcnologically, this discovery shares the 
following point with the subject of this paper: The non- 
classical decrease of the moment of inertia occurs even 
in the system in which the existence of the macroscopic 
Bose condensate is not expected. There is, however, the 
following difference. In a normal liquid helium 4 at the 
vicinity of T\, helium 4 atoms actually exist, and the size 
of their coherent wave function with zero momentum is a 
subject of the problem, whereas in a solid helium 4, the 
existence of a hypothetical moving boson, a zero-point 
vacancy, has not been firmly established, although it is 
normally assumed to exist. There are a number of prob- 
lems to be clarified in a solid helium 4. 



APPENDIX A: COMPARISON TO THE 
THERMAL FLUCTUATION 

Compared with the thermal-fluctuation theories, the 
formalism in Sec. 3 has the following differences. (1) In 
the fluctuation theory of a quantity x in question, the cor- 
relation function of deviations x — x from a mean value 
x, gives the susceptibility, an average of which is taken 
with the Gauss distribution. Owing to the flatness of 
the minimum of the thermodynamic potential near T c , 
the susceptibility has (1 — T/T c ) dependence. (Alterna- 
tively, the introduction of Green's functions representing 
thermal fluctuations leads to the similar (1 — T/T c ) de- 
pendence.) On the other hand, the susceptibility in Sec. 3 
is given as the correlation function of total quantities x 
as in Eq.(9). Its average is taken with the Bose- Einstein 
distribution. Hence, the result has no (1 — T/T c ) de- 
pendence. (2) Since thermal fluctuation is essentially a 
local phenomenon, one notes its influence already in the 
low-order terms of the perturbation expansion in which 
only a few particles participate. On the other hand, the 
change of the moment of inertia I z requires a transfor- 
mation extending to the whole liquid. Hence, the higher- 
order terms in which many particles participate plays a 
significant role. Only after the perturbation expansion is 
summed up to the infinite order as in Eq.(18), one no- 
tices the change of I z at a temperature in which Eq.(20) 
is satisfied. 

(3) Thermal fluctuations gradually increase as T — > T c , 
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finally diverging at T c . For the change of the moment of at T c , and I z smoothly agrees with the conventional value 
inertia, however, no physical relevant quantity diverges of 7° (1 — p s /p) as in Eq.(26). 
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